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The control of chemical reactions is a recurring 
theme in physics and chemistry. Traditionally, 
chemical reactions have been investigated by tun- 
ing thermodynamic parameters, such as temper- 
ature or pressure. More recently, physical meth- 
ods such as laser 1 or magnetic field 2 control have 
emerged to provide completely new experimental 
possibilities, in particular in the realm of cold col- 
lisions. The control of reaction pathways is also a 
critical component to implement molecular quan- 
tum information processing 3 . For these under- 
takings, single particles provide a clean and well- 
controlled experimental system. Here, we report 
on the experimental tuning of the exchange re- 
action rates of a single trapped ion with ultra- 
cold neutral atoms by exerting control over both 
their quantum states. We observe the influence 
of the hyperfine interaction on chemical reaction 
rates and branching ratios, and monitor the kine- 
matics of the reaction products. These investiga- 
tions advance chemistry with single trapped par- 
ticles towards achieving quantum-limited control 
of chemical reactions and indicate limits for buffer 
gas cooling of single ion clocks. 

The full control over all quantum mechanical degrees 
of freedom of a chemical reaction allows for the identifica- 
tion of fundamental interaction processes and for steer- 
ing chemical reactions. This task is often complicated 
in heteronuclear systems by a multitude of possible re- 
action channels, which make theoretical treatments very 
challenging. Therefore, focussing on the best-controlled 
experimental conditions, such as using state-selected sin- 
gle particles and low temperatures, is crucial for the in- 
vestigation of chemical processes at the most elementary 
level. The hybrid system of trapped atoms and ions offers 
key advantages in this undertaking. On the one hand, ion 
traps offer a large potential well depth in order to trap 
the reaction products for precision manipulation and in- 
vestigation. On the other hand, contrary to pure ionic 
systems, there is no Coulomb-barrier between the parti- 
cles which fundamentally prevents chemical reactions at 
low temperatures. Therefore, the efforts to control the 
motional degrees of freedom of one^~— and bother— reac- 
tants in hybrid atom-ion systems have paved new ways 
towards cold chemistry. The yet missing component is 
the simultaneous control of the internal degrees of free- 
dom. 

The interaction between an ion and a neutral atom at 
long distances is dominated by the attractive polarization 
interaction potential, which is of the form 



Here, C4 = a q 2 /(ineo) 2 is proportional to the neutral 
particle polarizability ao, q is the charge of the ion, eo is 
the vacuum permittivity, and r is the internuclear separa- 
tion. Inelastic collisions take place at short internuclear 
distances. In the cold, semiclassical regime this requires 
collision energies above the centrifugal barrier— ~—. Such 
processes are referred to as Langevin type collisions 
and happen, even for cold collisions 7,9 , at an energy- 
independent rate jLangevin = 2lTy/Ci/fl n Q . Here \i is 
the reduced mass of the collision partners and n a is the 
neutral atom density. More subtle effects, such as the hy- 
perfine interaction, which may lead to atom-ion Feshbach 
resonances, are not included in the polarization potential 
and have been investigated theoretically only so far—. 

Experimentally, reactive Langevin collisions in the po- 
larization potential have been investigated in ground 
state collisions 7 - " 1 1 1 2 1 1 8 , which have displayed relatively 
low rates for inelastic collisions, except for resonant 
charge exchange 7 .. Recently, first steps towards under- 
standing reactive collisions in excited electronic states 



have been made using large ion crystals 
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ing either a dominant contribution from very short-lived 
electronic states in the Rb+Ca + system-^ or, contrarily, 
a negligible contribution from excited state collisions in 
the Ca+Yb + system^. 

Here, we demonstrate how control over the internal 
electronic state of a single ion and the hyperfine state of 
neutral atoms can be employed to tune cold exchange re- 
action processes. We study quenching, charge-exchange, 
and branching ratios, and we use near-resonant laser light 
to control the rates. Our measurements show a large 
sensitivity of the charge-exchange reaction rates to the 
atomic hyperfine state, highlighting the influence of the 
nuclear spin on atom-ion collisions. 

In our experiment we study collisions between ul- 
tracold 87 Rb atoms and single 174 Yb + ions, for which 
1 Langevin /n a = 2.1 x 10 _15 m 3 /s. We start by determin- 
ing the inelastic collision loss rate coefficients for the long- 
lived 2 D 3 / 2 (radiative lifetime 52 ms) and 2 F 7 / 2 (radia- 
tive lifetime 10 years) states of the ion (see Figure 1 and 
Methods). The collisional stability of these states is of 
importance in buffer gas cooled ion clocks 19 and quantum 
information processing. By optical pumping, both states 
are prepared as "dark" states (see Methods), in order 
to study pure two-body collisions without the presence 
of light. This approach differs fundamentally from pre- 
vious experiments in atom-ior>ii excited state collisions, 
which always have been in presence of near-resonant laser 
light. We measure the inelastic loss rate 7^ by immers- 
ing the ion for a variable time t into the neutral atom 
cloud and determining the survival probability 8 ' 9 ' 20 ' 21 
P s = exp(— "f(t). In order to understand and model 
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FIG. 1. Level scheme and inelastic processes of the ion-neutral 
system, a) Left: Yb + level scheme with the transitions used 
for optically pumping the ion (to scale). Middle: Level scheme 
of a Rubidium atom. Right: Asymptotic level schemes of the 
two channels Yb + +Rb and Yb+Rb + . The collision is initi- 
ated in the Yb + +Rb manifold and the Yb+Rb + manifold can 
be populated by charge-exchange processes. Pictorial repre- 
sentation of the charge exchange reaction (b) and quenching 
from an excited state (c). Shown are the filled core electronic 
shells and the relevant valence electrons. 



our data, we make three assumptions: Firstly, the atom 
and the ion can undergo an inelastic process only in a 
Langevin-type collision. Secondly, the characteristic col- 
lision time is significantly shorter than the radiative life- 
time. Thirdly, inelastic collisions can only be exothermic 
thanks to the very low kinetic energy of the colliding part- 
ners and the absence of resonant laser light. We model 
It = e lLangevin as proportional to the Langevin collision 
rate using the state-dependent proportionality constant 
e. For the electronic ground-state of the ion interacting 
with the \F = 2,mp = 2) hyperfine ground state of Rb 
we find 9 f 



measure t% 2) = 1.0 ± 0.2 and ep 2) = 0.018 ± 0.004 (see 
Table HJ. The constant e = 1 corresponds to the largest 
allowed inelastic collision rate in the semiclassical model, 
and even in near-resonant charge- exchange between equal 
elements of atoms and ionsi it was expected and approx- 
imately found to be e = 1/2. In this regard, our results 
are unexpected since the 2 -D 3 / 2 state of Yb + combined 
with the 2 S , 1 / 2 -state of Rb is off-resonant by 0.14 eV to 
the next available asymptotic state in the Yb+Rb + man- 
ifold (see Fig. 1). Quite oppositely, the 2 -F 7 / 2 -state has 
nearby asymptotic states but its inelastic collision rate is 



significantly lower. 

We compare these results now to inelastic collisions in 
the absolute lowest hyperfme state \F — l,mp = 1) of 
the neutral atom. The hyperfine energy difference be- 
tween the \F — l,mF — 1) and the \F = 2,rriF = 2) 
states of Rb is 30 fJ.eV, which is larger than the collision 
energy but negligible on the scale of the molecular poten- 
tials or the trap depth (250 meV). We find a significantly 
enhanced probability for inelastic collisions with the ion 
in the electronic ground state e)^' 1 ^ = (35 ± 11) x e^' 2 \ 
which demonstrates an important role of the hyperfine 
interaction. For collisions with the ion in the -D3/2 state, 
the value of e remains unchanged. 

For collisions in the electronic ground state of the 
ion, exothermic inelastic collisions are inevitably asso- 
ciated with a chemical reaction, leading to charge ex- 
change or possibly the radiative association of molecules. 
The excess energy is converted into photons and/or ki- 
netic energy of the reaction products. Depending on the 
amount of kinetic energy released, the reaction products 
are ejected from the ion trap or kept and detected by 
in-trap mass spectrometry^ 2 ^. A partial conversion of 
internal energy into kinetic energy in the predominantly 
predicted-^ case of radiative charge exchange could be 
caused by the energy-dependent Franck-Condon overlap 
of the wave function of the entrance channels (A 1 S + and 
a 3 S+) and the exit channel of the (RbYb)+ poten- 

tials, owing to a 70% smaller polarizability of the X 1 E + 
state (assuming negligible momentum transfer by the 
photon and constant orbital angular momentum). We 
control the distribution of the kinetic energy release in 
the charge exchange reaction by employing the hyperfine 
state of the neutral atoms; for atoms in \F = 1, mp = 1) 
we retain Rb + with (48 ± 3)% probability in the ion 
trap and for atoms in \F = 2, mp = 2) we find it with 
(35±3)% probability. We also study inelastic collisions in 
the excited dark D— and F— states, where the probabili- 
ties for observing Rb + are quite similar to each other (see 
Table |T|). We do not observe the formation of (YbRb) + or 
RbJ molecules. The search for the latter was performed 
with the isotope 176 Yb + in order to achieve the required 
mass resolution. 

With regards to collisional quenching between differ- 
ent electronic levels of the ion, we are principally able to 
detect two different quenching scenarios: 2 F 7 / 2 — ► 2 £1/2 



= xlO 5±0,a , and for the excited states we 2 , 



and 2 £>3/2 — > 2 -FV/2- We directly observe quenching of 
"7/2 — ► 2 S1/2 as shown in Figure 2, where we exem- 



plarily show two different experimental runs, with and 
without quenching (see also Methods). After the quench- 
ing, we observe an increased temperature of the ions, 
which indicates a release of kinetic energy smaller than 
the depth of our trap, despite of the large energy gap be- 
tween the 2 Si/ 2 and the 2 -FV/2 states. The quenching rate 
from 2 i?3/2 — >• 2 i 7 V/2 was observed not to be detectable 
above our background rate. 

Having established the inelastic collision parameters of 
the metastable D- and instates without resonant laser 
light, we now turn our attention to inelastic collisions in 
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TABLE I. Measured proportionality constant e and branching ratios. The left part of the table refers to neutral atoms in the 
\F — 2, rap = 2) state and the right part to neutral atoms in the \F — 1, rnp — 1) state. Loss events of the charged particle are 
always chemical reactions for the S— states, whereas for D— and F— states quenching and chemical reactions can contribute. 
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FIG. 2. Collisional quenching from 2 i 7 V/2 to 2 Si/2- A two-ion 
Coulomb crystal is prepared in the 2 J7/2 dark state by optical 
pumping, then interacts with the neutral atoms, and is subse- 
quently probed by laser fluorescence (see Methods). An early 
appearance of laser fluorescence indicates collisional quench- 
ing (red curve) as compared to no quenching (blue curve). 
The large fluorescence dips indicate a high temperature of 
the ion crystal. 
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FIG. 3. Inelastic collision control by laser light, a State popu- 
lation measurement of the ion in the absence of neutral atoms 
as a function of the detuning of the repump laser at 935 nm. b 
Inelastic collision rate for the same detuning data in presence 
of neutral atoms at n a — 1 x 10 18 m -3 . The solid line shows 
the theoretical result of Eq. [2] c Density dependence of the 
inelastic collision rate with the repump laser on resonance. 
The exponent of the power-law fit (solid line) is 0.98 ± 0.02. 



the presence of laser light on both the 2 S 1 / 2 — 2 P\/ 2 (369 
nm) and the 2 D 3 / 2 — 3 D[3/2]i/ 2 (935 nm) transitions. 
The purpose of the light is to experimentally tune the 
rates and to observe the occurrence of inelastic collisions 
in real time. The radiative lifetime of the 2 P\/ 2 state 
(8 ns) is too short compared to the collision rate to pro- 
vide pure P-state measurements. Here and in the follow- 
ing, we assume that the instate is unoccupied since the 
411-nm light is turned off and collisional quenching into 
the instate has been measured to be negligible. There- 
fore, the overall inelastic collision rate is determined by 
a mixture of S-, P-, D-, and D[3/2]-states 

7£ = 2-KyJC±/ p,n a (p s e s +ppe P + p D e D + p D [ 3 / 2 ] e D[3 / 2] 

(2) 

Here, p x is the occupation probability of state x, which 
we determine experimentally for different settings of laser 
intensities and detunings (see Methods) . In Figure 3a we 
show the state populations of the S, P, D and 3 D[3/2]i/ 2 
states as we vary the frequency of the laser at 935 nm on 
the 2 D 3 / 2 — 3 D[3/2] 1 / 2 transition. Figure 3b shows the 
associated change in the reaction rates in the presence of 
the neutral atoms, which are proportional to the inelastic 
collision rates using the branching ratios of Table 1 . We 
demonstrate tuning of the reaction rate by one order of 
magnitude and we find it closely follows the model of 
equation (2) indicating a dominant contribution from the 
-D-state. From these data we also extract ep = 0.1 ± 
0.2, which is small and consistent with zero. Figure 3c 
shows the linear scaling of the inelastic collision rate with 
neutral atom density, confirming the picture of binary 
collisions. 

Finally, we turn our attention to the observation of 
the kinematics of the collision products. In Figure 4 we 
show the fluorescence at 369 nm during the interaction, 
with the 935 -nm light on. Subplots 4a-4c show individ- 
ual experimental runs, in which the initial sharp loss of 
fluorescence results from an inelastic collision. If some 
kinetic energy is released in this process, the ion is on 
a large trajectory in the trap (much larger than the size 
of the atom cloud) but can be re-cooled, which is sig- 
nalled by the relatively slow increase of the reoccurring 
fluorescence from which we determine a lower bound for 
the release energy of 8meV. We observe in 4% of our 
events that the fluorescence of the ion reoccurs after a 
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FIG. 4. Monitoring of inelastic atom-ion collisions, a-c 
Recorded fluorescence for selected events. After a quick initial 
loss, the fluorescence reoccurs after random times when the 
kinetic energy released in the inelastic collision has been re- 
moved by cooling. Then, a second collision occurs and the ion 
disappears again. The curves are vertically offset for clarity, 
d Sum of 343 repetitions of the experiment, which is fitted 
with an exponential decay for short times, e Zoom-in on the 
initial decay of plot d. The solid line is an exponential fit to 
the data. 



certain time, before it undergoes a second inelastic colli- 
sion and is dark again. As we keep the Yb + ion in the 
trap, these processes are not charge-exchange reactions 
but quenching processes with a kinetic energy release less 
than « 250 meV. Since the ion is colliding in highly ex- 
cited electronic states, such as P\ii and -D3/2 with inter- 



nal energies of « 3eV, this suggests a mostly radiative 
decay into the ground state Si/ 2 - We have ruled out that 
the reoccurrence events are linked to the dissociation of 
a potential molecular ion in a secondary collision with 
neutral atomsi^. To this end, we have performed this 
measurement also with extremely short interaction times 
between neutral atoms and the ion (on the order of a few 
collision times) and performed mass spectrometry on the 
dark ion after the collision confirming the absence of a 
molecular ion. 
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METHODS 
A. Preparation of ultracold atoms and ions 

We prepare 4 x 10 5 neutral 87 Rb atoms in the \F = 
2,mp — 2) hyperfine state of the electronic ground state 
at temperatures down to T sa 200 nK in a harmonic mag- 
netic trap of characteristic frequencies^ 3 - (ui x ,u) v ,oj z ) = 
2n x (8, 26, 27) Hz. By changing atom number and tem- 
perature, we tune the atomic density by two orders of 
magnitude. The atoms can be transferred into an op- 
tical dipole trap formed by two crossed laser beams at 
1064 nm. Here, the atoms can be transferred into the 
ground hyperfine state \F = l,mp = 1} using a reso- 
nant microwave pulse. At the same location, we trap 
single Yb + ions in a radio-frequency Paul trap with sec- 
ular trap frequencies of uj± — 2n x 150 kHz radially and 
LJax = 27r x 42 kHz axially^. We use standard techniques 
to cool and detect single ions or small ion crystals on the 
cycling transition S\/ 2 — 2 P\/2 near 369 nm wavelength 
(see Fig. 1). From the excited 2 P\/2 state, radiative 
decay populates the 2 D 3 / 2 state with a probability 2 ^ of 
i=3 1/200. For efficient laser cooling and detection, a laser 
at 935 nm pumps the population via the 3 -D[3/2] 1 / 2 state 
back into the cooling cycle. Preparation in 2 D 3 / 2 (radia- 
tive lifetime 52 ms) is achieved by optical pumping from 
the Si/2 state using laser light at 369 nm in less than 
10 fis in the absence of light at 935 nm, which is signif- 
icantly faster than the inverse of the Langevin collision 
rate of a few 10 3 s _1 . Alternatively, the 2 F 7 / 2 state (ra- 
diative lifetime 10 years) is populated by optical pump- 
ing on the 2 Si/2 — 2 D$/2 une with spontaneous decay 
on the 2 D 5 / 2 — 2 F- ! / 2 transition—. The ion is spatially 
overlapped with the center of the ultracold neutral atom 
cloud by displacing the two independent trapping poten- 
tials. The typical collision energies^, due to residual 
micromotion, are on the order of ~ 50 mK. 

B. Detection of collisional quenching 

We prepare two ions in a small Coulomb crystal at 
time t = in the 2 Sx/ 2 ground state, which we detect 
by near-resonant laser fluorescence at 369 nm (together 
with repumping at 935 nm) . The optical pumping into 
the 2 i r 7/2 state using light at 411 nm is switched on at 
t = 500 ms for 500 ms and typically within 100 ms the 
ion is pumped into the desired state (see Figure 2). As 
a result, the fluorescence counts on the 369-nm transi- 



tion drop to zero. After the interaction with the neutral 
atoms (effective duration typically 16 ms) and the subse- 
quent removal of the neutral atom cloud, the fluorescence 
on the 2 Si/ 2 — 2 P\/ 2 transition (together with repumping 
at 935 nm) is probed from t = 5000 ms onwards. If the 
ion is still in the 2 F 7 / 2 state it will not scatter photons, 
however, if it has been quenched to 2 5 , 1 / 2 , we observe flu- 
orescence. At t = 6000ms, the ions are optically pumped 
from the F-state back into the 5-state to ensure that no 
ions have been lost. The fluorescence count rate during 
the preparation part of the sequence is lower owing to the 
presence of an offset field from the magnetic atom trap, 
which is turned off when the atoms are released. 



C. Experimental determination of the ion 
electronic state occupation 

Owing to the experimental situation of a multilevel 
system in presence of a magnetic field and imperfect po- 
larization of the laser beams, we determine the state 
populations p x = t x /t c experimentally rather than re- 
lying on theoretical modelling. Here, r c is the average 
time between two spontaneous decays on the 2 P\/ 2 — 
2 D 3 / 2 transition, t x is the average time spent in state 
xe{S,P,D,D[3/2}} during r c , and r c = £ x r x . The 
time Tp is given by the decay rate T of the P\i 2 state 
and the branching ratio of the P-state into S— and 
D— states as Tp = 200/r. A sequence of alternating 
pulses of the 369 nm and the 935 nm lasers are used to 
determine t$ and Tp,. Tp + t$ is observed as the ex- 
ponential decay constant of fluorescence after the 369 
nm laser is pulsed on. The pulse length is set to get 
complete depletion of the S-P system into the D-state. 
For a given setting of intensity and detuning of the 935- 
nm laser, the average number of photons N 3 qq per 369- 
nm pulse depends on the duration of the preced- 
ing repump pulse. Varying the length £935 allows us to 
retrieve Tp and the photon detection efficiency 77 from 
the fit 200?y[l — exp(—t 935 / Tp)] = A 36 g to the counted 
photon number. We find 77 = (2.1 ± 0.1) x 10~ 3 , in 
accordance with the numerical aperture of the imaging 
system and the quantum efficiency of the single photon 
counter. The lifetime of the excited 3 D[3/2]!/2 state is 
r D[3/2] = 40 ns and the branching ratio is 98% into Si/ 2 
and 2% into D 3 / 2 . Decay into P\/ 2 is not dipole-allowcd 
due to parity 2 ^ and therefore we do not have to consider 
a cascade through intermediate levels. 



